Case: Schnitzer 5 

Serial No. 10/082870 A^f 
Filing Date: February 25, 2002 n 

< 12) INTERNATIONAL APPLICATION PUBLISH ED UNDER THE PATENT COOPERATION TREATY (PCT) 

(19) World Intellectual Property Organization 
International Bureau 

(43) International Publication Date 

20 June 2002 (20.06.2002) PCT 




(51) International Patent Classification 7 



(10) International Publication Number 

WO 02/48688 Al 



G01N 21/49 (72) 

(21) International Application Number: PCT/USOJ/48753 ^ ^ 



(22) International Filing Date: 

14 December 2001 (14. 12. 2001) 



(25) Filing Language: 

(26) Publication Language: 



English 
English (81) 



(30) Priority Data: 

60/256 J 19 



1 5 December 2000 ( J 5. 1 2.2000) US 



(63) Related by continuation (CON) or continuation-in-part 
(CIP) to earlier application: 

US 60/256,119 (CIP) 

F,,cd on 15 December 2000 (15.12.2000) 

(71) Applicant (for all designated States except US): THE 
GENERAL HOSPITAL CORPORATION [US/USj; 55 
Fruii Street, Boston, MA 02114 (US). 



Inventors; and 

Inventors/Applicants (for US only): BOAS, David, A. 

I US/USJ; Eight Maple Street, Newmarket, NH 03857 (US) 
DUNN, Andrew, K. jUS/US]; 10 Ridge Street, Arlington 
MA 02474 (US). ' 

(74) Agent: FASSE, Peter, J.; Fish & Richardson PC, 225 
Franklin Street, Bosion, MA 021 10-2804 (US). 

Designated States (national): AE, AG, AL. AM, AT AU 
AZ, BA, BB, BG, BR, BY, BZ, CA, Ci I, CN, CO, Cr! CU 
CZ, DE, DK. DM, DZ, HC, EE, ES, FI, GB, GD GE GH 
GM, HR, HU, ID, IL, IN, IS. JP, KE. KG, KP, KR, KZ, LC 
LK, LR, LS. LT, LU, LV. MA. MD, MG, MK, MN MW 
MX, MZ, NO, NZ, PL, PT, RO, RU, SD, SE. SG SI SK 
SL. TJ, TM, TR, TT, TZ, UA, UG, US, UZ VN Yu' ZA 1 
ZW. 

(84) Designated States (regional): ARIPO patent (GH GM 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW)' 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM ) 
European patent (AT, BE, CH, CY, DE, DK, ES, FI, FR 




oc 
oc 



o 




Detector/Pinhole 

' \ x- 



(57) Abstract: Indirect, three-dimensional images of features within mrhiH „t i 

p:« of radiative transport to provide higher resolution ZIZ Sn diS ST th ™° n ™ Cled based °" P*»cJ. 
enablc one to "see" into turbid or dense samples, such ™ tfss ue ,r lt"T 2 - T ^ melhods 

materials, to a depth of 50 pm to 10 mm or more, rttofoZ ^toZ^M *' "'T!,' """^ and 0(her 

source (1 1 ; illuminates a sample (18) and muitiplv W^^^'^^.'^T* P^' radiali °" tr0m 3 

one detector (22). "tterea optical radiation (20) emitted from the sample is detected by at least 



4 




r-R PR IP IT LU MC, NL. Y\\ Sli. TR), OAl'l palcnl 
J; S ep, CO Ci, CM. CjA, ON, «Q. OW, ML. MR. 
Nli. SN.TD.TG). 



Published: 

wiih international search report 



- before the expiration of the time limit for amending the 
claims and to be republished in the event of receipt of 
amendments 

For two-Utter codas and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begtn- 
ning of each regular issue of the PCT Gazette. 



WO 02/48688 



PCT/USO 1/48753 



INDIRECT MODE IMAGING 



CROSS-REFERENCE TO RELATED APPLICATION 

This application claims priority from U.S. Provisional Patent Application Serial No 
60/256,1 19, filed on December 15, 2000, which is incorporated herein by reference in its 
entirety. 

STATEMENT AS TO FEDERALLY SPONSORED RESEARCH 

This invention was made with Government support under NIH 1 P41 RR 14075 and 
NCI T32 CA09362 awarded by the National Institutes of Health. As a result, the 
Government has certain rights in the invention. 

TECHNICAL FIELD 

This invention relates to optical imaging, and more particularly to imaging of turbid 
or dense media. 

BACKGROUND 

Biomedical optical imaging has been studied over the past decade at both the 
microscopic and diffuse levels. Microscopic methods such as confocal reflectance (see e g 
Rajadhyaksha et al., J. Investigative Dermatol, 104:946, 1995) and multiphoton fluorescence 
microscopy (see, e.g., Denk et al., Science, 248:78, 1990) provide depth resolved, micron 
resolution images of tissue structure and function, but are limited to depths less than several 
hundred microns. Optical coherence tomography (OCT; see e.g., Huang et al, Science 
254:1 178, 1991) can penetrate to greater depths (about 1.0 mm), but at the expense of lower 
spatial resolution (5 - 20 urn). The limitations on the penetration depths of these methods 
anse from the fact that each relies on the detection of single scattered light, and since tissues 
are highly scattering, the probability of detecting single scattered light through several 
hundred microns becomes prohibitively small. 

Imaging with diffuse light on the other hand, typically yields spatial resolution on the 
order of several millimeters to a centimeter with penetration depths of several centimeters 
Since the detected tight has been scattered many times, an image reconstruction algorithm 
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m ust be employed to reconstruct scattering and absorption perturbations. Various 
.construction strategies have been used (see, e.g., Arridge et al., Inverse ProUens, 
1 5QYR41 1999) and many are based on the diffusion approximation to the Boltzmann 



transport equation. < 

The spatial regime lying between that covered by the microscop.c and drfS.se 
techniques remains relatively unexpired. This is in part dueto.be complex of the 
description of light interaction with tissue on a length scale of only a few soattermg lengms. 
While microscopic methods yield direct reflectance or fluorescence maps of the tissu 
architectiue, diift.se methods utilize a series of indirect measurements and subsequently 
re— an image by solving the inverse problem In this intetmedtate spanal regrme, 
direct imaging of single-scattered light is not feasible. 

SUMMARY 

The invention is based on the discovery that induect. three-dimensional images of 
features within turbid or dense samples, such as tissues in living animals or humans, 

radiative nansport to provide higher resohation images than diffuse .magntg mefltoda, whale 
enabling penetration depths significantly greater ton microscopic imaging methods. 

to general, the invention features methods of indireo, mode imaging of a sample by 
(a) ifluminating the sample with optical radiation from a source; (b) receiving optical 
2 „ radiation emitted from fire sample with one or more detectors-, (c) dtgitizing tire optica, 
radiation received fiomthe detectors to generate a digitized signal and transmuting the 
digitized sigrtal to a processor; and (d) processing the digitized signal to reconstruct an mtage 
of spatially varying optica! properties of feamres within tite sample by performing a non- 
linear minimization between the digitized data and a transport-based photon nugraoon 

m0de '' m these methods, a "non-linear minimization" is an algorithm designed to minimize 
,he difference between Wo quantities, as described in fiuther detail below. A Wport- 
hased photon migration model" is a model that describes the propagation of hgh, (photons) 
toough a turbid medium such as biological tissue, as described in further detail herem. 
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In these methods, the detectors can each be located at a position offset from the 
source, and one or more of the offsets of each detector-source pair can be different. The 
image S Ma (r) can be reconstructed by calculating a matrix equation of the form y = Ax, 
wherein Ma is the absorption coefficient, r is a position within the sample, x = 5 Ma , y is the 
vector of detector signals for each source-detector pair, and^, is an element of matrix A 
representing the product of the Green function for the transport equation (G) and the first- 
order (background) detector signal (Lo) for each measurement / at each voxel j. In these 
methods, tire optical radiation can be visible, near-infrared, or other radiation, and there can 
be one or more sources. Alternatively, the optical radiation can be scanned from the source 
across the sample to simulate multiple sources. For example, 10, 15, 20, 30, 40, 50 or more 
detectors (or sources) can be used. Tire offset can be, for example, from 0.1 to 10 mm. 

In another aspect, the invention features systems for indirect imaging of a sample that 
includes (a) a probe comprising a source optic fiber, one or more detector optic fibers, and a 
distal end, wherein a distal end of the source optic fiber and each of the detector optic fibers 
extends through and ends in the distal end of the probe; (b) an optical radiation source 
connected with a proximal end of the source optic fiber; (c) one or more photodetectors, each 
connected to a proximal end of one of the detector optic fibers, to detect and convert optical 
radratron from each detector optic fiber into a digital signal corresponding to the light emitted 
from the sample; and (d) a processor that processes the digital signal produced by the 
Photodetectors to provide on an output device an image of spatially varying optical properties 
of features within tire sample, wherein the image is reconstructed by performing a non-linear 
mrnrmrzation between the digitized data and a transport-based photon migration model. 

In these systems, the distal end of each detector optic fiber can be offset from the 
drstal end ofthe source optic fiber, e.g., by about 0.1 to 10mm. In addition, in this system 
the processor is programmed to process the digital signal to provide an image S M r) that is 
reconstructed by calculating a matrix equation ofthe form y ^ Ax, wherein » 3 is the 
absorption coefficient, r is a position within the sample,. = S Ma , y is the vector of detector 
srgnals for each source-detector pair, and ,4, is an element of matrix A representing the 
product ofthe Green function for the transport equation (G) and the first-order (background) 
detector signal (Lo) for each measurement i at each voxel / 
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As above, the one or more of the offsets of each detector-source pair can be different, 
the optical radiation can be near-infrared radiation, and tire optical radiation can be scanned 
from the source optic fiber across the sample to simulate multiple sources. In some 
embodiments, 10, 15, 20, 30, 40, 50 or more detectors (or sources) can be used. The offset 
can be from 0.1 to 10 mm, and can vary form one to the other offset. 

In another aspect, the invention features a probe for indirect mode imaging. The 
probe includes a source optic fiber, one or more detector optic fibers, and a distal end, 
wherein a distal end of the source optic fiber and each of the detector optic fibers extends 
through and ends in the distal end of the probe, and wherein the distal end of the source optic 
fiber is offset from each distal end of the detector optic fibers by 0.1 to 10 mm. The probe 
can further include a scanning mirror to scan optical radiation across a sample. The probe 
can have a variety of different length offsets for the one or more source-detector pairs. 

Unless otherwise defined, all technical and scientific terms used herein have the same 
meaning as commonly understood by one of ordinary skill in the art to which this invention 
belongs Although methods and materials similar or equivalent to those described herein can 
be used in the practice or testing of the present invention, suitable methods and materials are 
described below. All publications, patent applications, patents, and other references 
mentioned herein are incorporated by reference in their entirety. In case of conflict, the 
present specification, including definitions, will control. In addition, the materials, methods, 
and examples are illustrative only and not intended to be limiting. 

The new imaging methods provide the ability to "see" into a sample, e.g., tissue, to a 
depth of 50 urn to 10 mm, which is far greater than is possible using microscopic imaging 
methods, which typically provide a view of only the surface of a sample. At the same time, 
the new methods provide a resolution of 10 urn to about 1 mm, which is significantly greater 
than is possible using diffuse imaging methods that provide useful penetration depths, but a 
resolution of only 5 to 10 mm. The new methods can be used in various diverse applications 
such as biomedical applications including functional imaging of the cerebral cortex through 
an intact skull, endoscopic imaging of small lesions too deep for conventional microscopic 
techniques to image and too small for diffuse methods to resolve, and ophthalmologic^ 
, imaging, e.g., of layers within the retina. Other applications include imaging of ceramics, 
semiconductors, and other turbid or dense materials, e.g., during processing. 
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A S1 gnificant advantage of the technique described herein is that it enables for the first 
time optical imaging in the intermediate regime between microscopy and diffuse optical 
tomography. Other features and advantages of the invention will be apparent from the 
following detailed description, and from the claims 

BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a graph illustrating the trade-off between penetration depth and resolution 
for vanous imaging methods including the new indirect mode imaging methods. 

FIG. 2 is a schematic diagram of a general imaging geometry useful in the new 
imaging methods. 

FIG. 3 is a graph of a perturbed signal computed with a first Bom approximation 
(solid lines) and with a full perturbative Monte Carlo simulation (symbols) at source-detector 
separations of 500 um and 2.0 mm for an object located 1 .0 mm beneath the surface with Su. 
= 0.1 mm* 1 . 

FIG. 4A is an end view of an illumination/detection probe useful in the new imaging 
methods. e 

FIG. 4B is a schematic diagram of an illumination/detection probe in a scanning 

system. 

FIG. 5A is a flowchart illustrating the steps used to acquire and process optical signal 
data in the new methods. 

FIG. 5B is a flowchart illustrating the steps used to process the optical signal to 
reconstruct an image of a feature in the sample. 

FIG. 6 is an endoscopic illumination/detection probe for use in the new methods. 

FIG. 7 is a computer-based simulation of the distribution of photons for a source- 
detector separation of 750 um. 

FIG. 8 is a graph of normalized radial intensity at z = 500 um vs. position on the x- 
axxs for the cases where the angular dependence is considered ("transport" - solid line) and 
the case where it is ignored ("diffusion" - dashed line). 

FIGs. 9A and 9B are a pair of reconstructed images of a 100 um absorbing object 
located at a depth of 1 .0 mm (9A) or 2.0 mm (9B) below the surface using simulated data 
(Su a = 0.1 mm" 1 ). 
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reference symbols in the vanous drawings indicate like elements. 



DETAILED DESCRIPTION 

The invention is teed on the discover <ha, indirect, taee-dimensiona! images of 

uansport to provide higher resolution images than difmae tmagtng mefcods, 

penelion depths aigniucanrly greater titan microscopic imaging methods. The new tndtrec, 

hving animais orhumans, ceramics, plastics, hqnids, and other matenaU, to a depth of *ou. 

, resolutton man diffnse imaging methods, which provide somewhat greater penetranon 
depths, but low resolution. 

r^npral Methodology . 
, 5 hot there is a significant gap in me ability to image features within a sample a, a depth 

There has always heen a trade-off between penenation depm and resolntton, hut to date, mere 

Thts uade-off is summarized in the graph in FIG. 1 . Althongh microscopy ^^iZ 
„ resolution (1 to .0 «m), it cannot he used to image much helow me surface of a sample (300 
Oumatbest) On the other hand, drfiuse mtagmg methods based on photon mtgratton 

2 about 1 cm. Optica, Coherence Tomography (OCT) has enabled ~ 
depths than micmscopy (up to about 0.5 to 1.0 mm), bu, still doea no, provtde an tmage o 
26 deeper features in a sample, and at the maxhnum penetratton depth provtdes a resolutton of 

^ 'iTew indirect mode rmaging 0MQ methods ate based on me principle that images 
can be reconstructed baaed on the frameworic of the radiative transport elation OttE) (see, 

e.g„ fshimaru e, a,., Wa~ ***** - — g in W» ™™"££L 
30 Press inc. 1978). An important distinction between reconstruction based on the R.TE and the 
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photon diffusion equation (DE) is the treatment of the angular dependence of the scattered 
light within the tissue. The new methods include both indirect microscopic imaging and 
indirect macroscopic imaging. 

FIG. 2 illustrates a general geometry 10 used to image in the new IMI methods, in 
which light is multiply scattered, but not yet diffuse. A source 1 1 directs light to a mirror 12 
and through lens 16 into target 18. Light reflected from the sample along the optical axis 
(beams 14) hits point 13. By translating a pinhole 21 and a detector 22, which is typically 
used in confocal reflectance measurements, to positions laterally displaced from the optical 
axis (at point 13), the pinhole 21 is imaged onto the detector 22 to a point laterally displaced 
from the source. Therefore, the amount of lateral displacement of the detector aperture from 
the source aperture determines the effective source-detector separation. This separation 
distance can vary from 0 to about 10 mm. Because the source and detector no longer lie in 
conjugate image planes, the detected light has been scattered at least one time. This scattered 
light is represented by dotted lines 20 in FIG. 2. 

As the source-detector separation distance increases, the detected light will probe a 
deeper and larger region of the sample. To maximize the sensitivity of the measurement to a 
particular region of the sample, several parameters can be varied in the geometry of FIG. 2. 
For example, the source-detector separation (aperture offset, 0 to about 10 mm, e.g., 1, 2, 3, 
5, or 7 mm), numerical aperture (0 to about 1.0, e.g., 0.1, 0.2, or 0.3), depth of focus intothe 
medium (0 to about 5 mm, e.g., 1, 2, 3, or 4 mm), wavelength (typically about 400 to about 
1500 nm, e.g., 570 to 650 nm or 1300 to 1500 nm), and pinhole diameter (about 100 um to 
about 1 .0 mm, e.g., 100, 150, 200, 250, or 300 urn, or 0.5, 0.75 or 0.9 mm can all be varied). 
In general, the size of the pinhole aperture will be larger than that used in confocal 
reflectance measurements since it will be more important to maximize the signal intensity 
than to provide a large degree of spatial filtering. 

Direct imaging of single-scattered light using the geometry of FIG. 2 is not feasible. 
Therefore, to form an image, an inverse problem must be solved as in the diffuse regime. 
Since the diffusion approximation is no longer valid due to the close proximity of the source 
and detector, a description based on the transport equation (see, e.g., Ishimaru et al., Wave 
Propagation and Scattering in Random Media, Vol. 1 (Academic Press, Inc. 1978)) is used. 
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To reconstruct an image of an absorption perturbation with a source-detector 
separation of only a few scattering lengths, one can use a linear reconstruction aigortthm 
based on the firs, Bom applanation .o the radiative ttansport elation. In the first Bom 
approximation, the absorption coefficient « and scattering coefficient W are ™tten as a 
6 sum of a homogeneous background component and a spatially varying pertuibatton, 

(r) = + *'« and ^ (r) " ^ + * ,(r) 

The radiance (i.e„ the amount of fight received at the detector), L, is then expressed as 
,„ , n esumofabackgroundandapermrbedcom P onen t Z, = L, + L,. For an absorptton 



15 



perturbation the zeroth- and first-order terms, U and L\, are given by: 
L a (r„ ,0. d ) = fj5(r, Q)G(r, Q | r„ , Q„ VKid'r 
i,(r. A) — fJ*.(r)I.(r..r,ft)GOr. r rf A)<^ 



(1) 
(2a) 



where G(r,Q|r rf A) is me Green function solution for the transport equa^on for a 
particular detector configuration. Perturbations in the scattering coefficient and phase 
function can be computed using the same approach and the first-order term for a scattenng 
20 perturbation is given by 

X l(rd ,a rf ) = -n^(rfc(r,r,Q)- fl.(r..r.^aQ^^|r, <W 

Eqs 2a and 2b illustrate that one of the primary differences between the new indirect 
25 mode imaging methods and imaging methods based on the diffusion approximation hes in 
the angular dependence of the radiance. Here, we only consider reconstruct of an 
absorbing object using Eq. 2a, but the same procedure can be followed to reconstruct a 
scattering perturbation. To reconstruct scattering perturbations, one uses exactly the same 
procedure but substitutes Eq. 2b for Eq. 2a. To simultaneously reconstruct an image of both 
30 absorption and scattering, one uses the summation of Eq. 2a and Eq. 2b. 



-8- 



WO 02/48688 

PCT/US01/48753 



10 



15 



20 



25 



30 



To reconstruct an image of & a (r) from Eq. 2a, L 0 and G are computed for a 
homogeneous background of known optical properties using a Monte Carlo simulation in a 
focused beam geometry (Dunn et al., Applied Optics, 35:3441 (1996). The Green function is 
computed by simulating the propagation of light from the detector to the medium and then 
utilizing reciprocity to determine the fraction of light reaching the detector from direction Q. 
at point r in the medium. Once the background radiance and Green function have been 
computed for all source-detector pairs, an image of S^{r) is reconstructed using die 
measured values, L u which are obtained here for simulated data. 

Before using Eq. 2a to reconstruct an image, the validity of the first Bom 
approximation for the transport equation must be established. To test the validity, a perturbed 
signal computed using Eq. 2a was compared with a perturbed signal computed using a 
perturbative Monte Carlo simulation of an absorbing object embedded in a homogeneous 
background (Sassarolia et al., Applied Optics, 37:7392, 1998). The absorbing object used in 
the simulated comparison was a cubic object with dimensions of 1 00 urn located at a depth 
of 1 mm beneath the surface of a turbid sample. The optical properties programmed for the 
simulation of the background and perturbation were „, = 0.001 mm-', Ms = 10 mm- 1 , g = 0.9, 
and S M * = 0. 1 mnf' . The perturbed signal computed with Eq. 2a and with the full ' 
perturbative Monte Carlo model as the absorbing object was laterally translated through the 
sample areplotted in FIG. 3, which shows the comparison of the perturbed signal computed 
with the first Bom approximation (solid lines) and with a full perturbative Monte Carlo 
simulation (symbols) at source-detector separations of 500 urn and 2 mm for an object 
located 1 mm beneath the surface with S Ma = 0.1 mm" 1 . A decrease in the perturbed signal is 
observed while either the source or detector is scanned across the object position. Although a 
useful image can be reconstructed using one source-detector pair at a fixed separation 
distance (offset), image resolution improves as additional source-detector pairs, at different 
offsets, are added to the calculation. 

The amplitudes of the perturbed signals in FIG. 3 have not been normalized and 
demonstrate that the absolute amplitudes of the perturbed signals are accurately predicted 
using the linear perturbation model. 

The comparison ulustrates that the first Born approximation accurately predicts the 
perturbed signal and that Eq. 2a can be used to reconstruct an image of « r) . Since the first 
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Bom approximation is a linear approximation, the magnitude of the perturbed signal is 
directly proportional to the magnitude of Based on the simulations, the linear 
approximation begins to deviate from the value predicted by the full perturbative Monte 
Carlo model at a of about 1.5 mm" 1 for a 100 urn object located at a depth of 1 mm. In 
5 general, the value at winch the linear approximation breaks down will vary depending on the 
size and depth of the perturbation. Results indicate that the maximum *. increases as the 
depth of the object increases and decreases as the physical size of the perturbation increases, 
in agreement with standard perturbation theory. 

10 Indirect Mod * Tmagine Devices 

The new indirect mode imaging systems include (i) an illumination and detection 
device such as a standard scanning confocal microscope (e.g., a Zeiss, LSM 410®), or an 
endoscopic probe, with various components that are used to illuminate a sample and collect 
light emitted from the sample, such as from tissue in a patient or animal body; (ii) a light 
1 5 source that is connected to the illumination and detection device by optic fibers; and (iii) an 
apparatus that converts the light into a digital signal and includes a processor programmed 
with algorithms that can process the digital signal into indirect, three-dimensional images 
that provide information about features within a sample, such as diagnostic and prognostic 
information about a tissue sample in a living animal or human patient. 
20 The new methods can be carried out with devices that are similar in some respects to 

confocal microscopes, but with a variable offset between the light source and one or more 
detectors. For example, the system can include an apparatus that is connected to a standard 
scanning confocal microscope. The light source can be a laser, such as a diode laser. Other 
light sources include light emitting diodes and lamps. Both visible and near-infrared light 
25 can be used, e.g., in the wavelength range of 400 to 1500 nm, but the choice of wavelength 
depends to some extent on the nature of the sample and the features to be imaged within the 
sample For example, to image hemoglobin within blood vessels, the preferred wavelength 
range is about 570 to 650 nm. To image ceramics, wavelengths of from 1300 to 1500 nm can 
be used to nunimize scattering from the sample. The light can be from a continuous wave 
30 (CW) source with a set amplitude and frequency, or can be modulated in either or both 
amplitude and/or frequency. In addition, the light can be pulsed light. Incoherent light 
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sources such as mode-locked solid-state pulsed lasers, and super-luminescent diodes can 
provide higher spatial resolution by providing optical measurements that are path-length 
resolved. 

FIG. 4A is a schematic end view of a confocal mumination/detection arm or probe 30 
useful in the new imaging systems. Probe 30 utilizes a modified form of a traditional 
confocal microscope. Probe 30 is connected to a light source 32 such as a laser, and multiple 
photodetectors 34, such as avalanche photodiodes (APDs) or photomultiplier tubes (PMTs) 
with different offsets from the optical axis. This configuration can be achieved through an 
array of optical fibers 36, 38 collected in a distal endplate 39. Source fiber 36, which can be 
located in the center of the probe or can be arranged to one side, delivers the light to the 
sample and also serves as a traditional confocal pinhole. The remaining detector fibers 38 
serve as off-axis pinholes and each is coupled to a photodetector 34. The fiber diameter and 
amount of offset can be varied to select signal level and penetration depth within the 
medium. For example, the average penetration depth of the detected light is about 0.5 mm 
for a source-detector separation of 750 mm and a fiber diameter of 100 mm. In addition, the 
offsets for some detector fibers to the source fiber can differ from the offsets for other 
detector fibers to the source fiber. In fact, the best images are obtained by having a 
multiplicity of source-detector separations. 

The new instruments 30 can be used in an imaging system 40 shown in FIG. 4B. As 
shown in FIG. 4B, Ulumination/detection probe 30 directs light to mirror 44 through lens 42, 
and then from the mirror through a second lens 46. The light continues through x/y scanning 
mirrors 48 and lens 50 and is scanned in a user-specified pattern (typically a raster scan) in 
the back focal plane of the objective lens 54. This scanning process is controlled by 
processor 56 using standard scanning confocal microscope scanning protocols. A portion of 
the light passes through a beam-splitter 52 and to objective 54, which focuses the light on to 
sample 18. Light reflected from the sample, returns along the same path, but some of the 
reflected light is directed to an imaging device 55, such as a charge-coupled device (CCD) or 
camera, that provides an image of the surface of the sample. The direct image of the surface 
of the sample can be co-registered or superimposed over the image provided by the new IMI. 
The remainder of the reflected light passes back to the illumination/detection probe 30, to be 
processed to reconstruct an image of features within sample 18 by processor 56. 
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Probe 30 of FIG. 4A is used with an apparatus that includes a standard analog-to- 
digital converter and a processor, as described in further detail below. Both the A-to-D 
converter and the processor can also be in a separate PC. As shown in FIG. 4B, such a 
processor 56 generally includes an input/control device 57, a memory 58, and an output 
device 59. The processor can be an electronic circuit comprising one or more components. 
The processor can be implemented in digital circuitry, analog circuitry, or both, it can be 
implemented in software, or may be an integrated state machine, or a hybrid thereof. 
Input/control device 57 can be a keyboard or other conventional device, and the output 
device 59 can be a cathode ray tube (CRT), other video display, printer, or other image 
display system. Memory 58 can be electronic (e.g., solid state), magnetic, or optical. The 
memory can be stored on an optical disk (e.g., a CD), an electromagnetic hard or floppy disk, 

or a combination thereof. 

As shown in the flowchart of FIG. 5A, the processor controls data acquisition and 
processing. First, in step 100, the processor controls the scanning of the mirror 44 to 
illuminate the sample, then, optionally, amplifies the detector signals to an appropriate level 
with analog circuitry. In step 1 10, the analog signal is digitized with an analog-to-digital 
converter, which may be contained within a computer. The digitized detector signal is stored 
in computer memory or on a hard disk. In step 120, the digitized detector signals are 
processed using an image reconstruction software package, which is described in further 
detail below. The resulting three-dimensional image is then displayed, e.g., in two- 
dimensional form (step 130). 

To reconstruct an image from the digitized detector signals, Eqs. 2a and/or 2b must be 
solved for the terms Sp a and/or 5* at each voxel location. The flowchart in FIG. 5B 
illustrates the steps in data processing step 120. The digitized data is processed by 
expressing Eq. 2a as a matrix equation of the formy =Ax, wherex = is the vector of 
optical properties at each voxel, y is the vector of detector signals for each source-detector 
pair, and Ay is an element of matrix A representing the product of the Green function for the 
transport equation (G) and the first-order (background) detector signal (Lo) for each 
measurement i at each voxel j. The Green function is obtained through Monte Carlo 
simulations of the light propagation in the medium for a particular source-detector 
configuration (Dunn et al., Applied Optics, 35:3441 (1996)). This matrix equation can then 
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be solved using standard techniques (see Arridge Inverse Problems, 15{2}:R41, 1999), such 
as truncated singular value decomposition. 

In step 121, the processor makes an initial "guess" of the optical properties (ju Sl Ma , g) 
for every image voxel. This list of optical properties at each voxel is definedas vector x. 
Next, in step 122 the processor calculates a predicted detector signal measurement, L p (x) 
given x. In step 123, the processor queries whether residual L p -L is above or below a 
threshold level typically set at the measurement noise level of the system. The measurement 
noise should typically be less than about 0.1 % of the measurement signal (L p -L/L < 10" 3 ). 
If the residual is less than the threshold, then the image is displayed (step 130). If the 
residual is larger than the threshold, then the mauix A is calculated given x (step 124). In 
step 125, the processor calculates an update tox given A (using any of several known 
techniques) and the difference between the measurement y and predicted measurement y p and 
uses the updated* to calculate another predicted detector signal measurement, L p (x) given 
the new x (step 122). 

The process steps 122, 123, 124, and 125 are repeated until the residual L p - L is 
small enough, in which case an image is displayed (step 130). 

FIG. 6 shows an endoscopic illumination/detection probe 60 for use in the new 
methods. In this probe 60, a laser 62 directs light through a beam-splitter 64, and through a 
pair of scanning mirrors 65. Light is focused through lens 65a onto the face of a fiber optic 
imaging bundle 66. A focused spot of light is scanned across the fiber bundle to couple light 
into a single fiber at a time. Fiber bundle 66 can be inserted into an endoscope 67 and a 
gradient index (GRIN) lens 68 can be used to focus the fight exiting each fiber within the 
bundle to a point inside a sample, e.g., inside tissue. Other combinations of lenses, such as 
micro-lenses and holographic lenses, can be used to focus the light within the tissue. 

Reflected light is also collected by GRIN lens 68 and is coupled back into the fiber 
bundle 66. The light exiting the bundle is imaged onto detector array 69 via lenses 65a and 
64a, and through scanning mirrors 65 and beam-splitter 64a. The detector array 69 consists 
of one or more detector elements whose position is laterally displaced from the optical axis. 
The amount of displacement determines the penetration depth of the detected light. For 
example, a source-detector separation of 750 urn provides a penetration depth of about 500 
um. This displacement can be varied to change the penetration depth of the light as it passes 
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m ,„ ,he tissue. By increasing the separation distance, the penetration depth increases (bn, tire 
resolution decreases). 

A pplications 

The new methods and devices provide the ability to image through several 
millimeters of turbid samples, such as human and animal tissue, with a resolution of a few 
hundred microns, and can therefore be used to address many new biomedical problems and 
applications. 

For example, the new methods and devices can be used to image subsurface blood 
vessels e g through the skin or tissue, located several millimeters within tissues since the 
size of the vessels is a few hundred microns and the difference in absorption between the 
blood vessels and the surrounding tissue will be about 1 mm-' in the near infrared. In another 
example, the endoscopic probe can be used to image walls of larger blood vessels, as well as 
the tissue outside the blood vessel walls. 

Another use of the new methods is in optimizing the use of the indirect mode of the 

scannmglaseroph^^ v 
sub-retinal structures where an annular aperture is used in place of a confocal aperture so that 
multiply scattered light is detected. Currently the use of the indirect mode is based on 
observation rather than a theory or model based on light scattering in tissue. Therefore, the 
>o model presented in this paper could be applied to optimize its use, because it provides a 

method for quantifying spatially varying tissue optical properties (absorption and scattering). 

Quantitative knowledge of these optical properties may reveal useful diagnostic informahon 

about the tissue state. . 

25 Software Im plementation 

The processing of the digital data corresponding to the light emitted from the sample 
can be implemented in hardware or software, or a combination of both. The method can be 
implemented in computer programs using standard programming techniques following the 
methods, equations, and figures described herein. Program code is applied to enter data to 
30 perform the functions described herein and to generate output information. The output 
information is applied to one or more output devices such as a display monitor. 
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Each program is preferably implemented in a high level procedural or object oriented 
programming language to communicate with a computer system. However, the programs 
can be implemented in assembly or machine language, if desired. In any case, the language 
can be a compiled or interpreted language. 

Each such computer program is preferably stored on a storage medium or device 
(eg., ROM or magnetic diskette) readable by a general or special purpose programmable 
computer, for configuring and operating the computer when the storage media or device is 
read by the computer to perform the procedures described herein. The computer program 
can also reside in cache or main memoiy during program execution. The processing methods 
can also be implemented as a computer-readable storage medium, configured with a 
computer program, where the storage medium so configured causes a computer to operate in 
a specific and predefined manner to perform the functions described herein. 

EXAMPLES 

The invention is further described in the following examples, which do not limit the 
scope of the invention described in the claims. 

Example 1 - Study of Angular Dependence of Radiance 

To study the angular dependence of the radiance, the integrand of Eq. 12 above was 
computed using a Monte Carlo simulation. The geometry for the simulation consisted of a 
beam focused at a depth of 500 urn beneath the surface with a numerical aperture of 0.4. The 
distribution of photons within the sample was computed for a sample with the following 
optical properties: ^ = 0.01 mm".*- 10*^ = 0.9, andn = 1.0. The source-detector 
separation was set at 750 um. 

The Green function G(r, Q | r d ,Cl d ) was computed from the results of the simulation 
by translating the computed photon distribution radially and utilizing the reciprocity of the 
photon paths, G(r,6\ r'M') = G(r',~Q | r ,Q) . 

As an example, FIG. 7 shows the photon distribution within the medium for photons 
reaching the detector, 



15- 



BNSDOCID <WO 02486HRA1 I > 



\N O 02/4X688 



PCT/US01/48753 



20 



jG(r„njr,fi)G(r,fl|r 1 „n (( )dC2 



(3) 



By examining the radial distribution of photons within the medium, the angular 
dependence can be illustrated as in FIG. 8, in which the radial distribution is plotted at a 
depth of 500 um for the angular resolved case (Eq. 3; solid line labeled "transport"), and the 
' when the angular dependence is ignored (Eq. 4; dash-dotted line labeled "diffusion"), 



case 



JG(7- , Q s | r, 0)dQ f G(r, Q | r d , Cl d )dQ 



(4) 



,0 The intensity in FIG. 8 has been normalized to the peak value in each case for 

comparison, but the absolute intensity is approximately an order of magnitude smaller for the 
case where the angular dependence is accounted for (Eq. 3) than where the angular 
dependence is ignored (Eq. 4). As shown in FIG. 8, the dip in intensity between the source 
and detector locations is more pronounced when the angular dependence is considered. 

1 5 Therefore, when reconstructing an image of an absorption inhomogeneity, the angular 
dependence of the photon distribution should be taken into account. 

Example 2 - Simulation of Image R econstruction 

Once the validity of the first Bom approximation has been established (FIG. 3), «i 
image of «r) can be constructed by writing Eq. 2 as a matrix equation of the form y = Ax, 
where y is the set of measurements and* - «r). Truncated singular value decomposition 
was used to invert this set of equations and find «r) using simulated data (Arridge, Inverse 

Problems, 15(2):R41, 1999). 

Reconstructed images of a 100 urn absorbing object = 1 mm- 1 ) at depths of 1 and 
9 nan are shown in FIGs 9A and B. The background optical properties of the medium were 
the same as those described in conjunction with FIG. 3 above. The set of measurements used 
in the simulated reconstruction consisted of source-detector separations ranging from 400 um 
to 2 mm in 200 um increments at numerical apertures of 0.2 and 0.4, for a total of 18 source- 
detector pairs (9 with a source and detector NA of 0.2 and 9 with an NA of 0.4). Each pair 
was translated 1.25 mm across the surface of the sample in 51 steps of 25 um yielding a total 
of 91S measurements. The focal point of the. source and detector was set to 1 mm beneath 
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the surface for all measurements. The singular value spectrum was truncated at 250 and this 
number was determined by considering a potential signal to noise ratio of 1 0 3 for the 
measurements (although in the present simulation there is no noise). The system was 
assumed to be shot noise limited and the number of singular values was chosen such that the 
magnitude of the perturbed signal was greater than the measurement noise in the total 
detected signal (background + perturbation). The images clearly indicate that lOOumaxial 
and lateral resolution is maintained to a depth of 1 mm (FIG. 9A) and 2 mm (FIG. 9B). 

Based on the images of FIGs 9A and 9B, it is clear that this method is capable of 
imaging in the spatial regime lying between the microscopic and diffuse regimes. 

OTHER EMBODIMENTS 

It is to be understood that while the invention has been described in conjunction with 
the detailed description thereof, the foregoing description is intended to illustrate and not 
limit the scope of the invention, which is defined by the scope of the appended claims For 
example, the principles described herein can be applied for imaging other types of radiation 
that has been multiply scattered and detected, i.e., imaging at other wavelengths such as x- 
rays and microwaves given proper illuitunation/detection devices. Other aspects, advantages, 
and modifications are within the scope of the following claims. 
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WHAT IS CLAIMED IS: 

I! A method of indirect mode imaging a sample, the method comprising 

(a) illuminating the sample with optical radiation from a source; 

(b) receiving optical radiation emitted from the sample with one or more detectors; 

(c) digitizing the optical radiation received from the detectors to generate a digitized 
signal and transmitting the digitized signal to a processor; and 

(d) processing the digitized signal to reconstruct an image of spatially varying optical 
properties of one or more features widnn the sample by performing a non-linear 
minimization between the digitized data and a transport-based photon migration model. 

method of claim 1 , wherein the detectors are each located at a position offset 



2. The 
from the source. 

3. The method of claim 2, wherein one or more of the offsets of each detector-source 
1 5 pair are different. 

4. The method of claim 1, wherein the image ^ a (r) is reconstructed by calculating a 
matrix equation of the formy = Ax, wherein 

fi z is the absorption coefficient, 
20 r is a position within the sample, 

x = Sfi a , 

y is the vector of detector signals for each source-detector pair, and 
A 9 is an element of matrix A representing the product of (i) a Green function: 
G(r , Q | r d , fi rf ) , for a transport equation (G): ji* (r) - A + W') ^ = ^ + 
25 ^ s (r), and (ii) a first-order, background detector signal (U), for each measurement i at each 
voxel j. 

5. The method of claim 1, wherein the optical radiation is near-infrared radiation. 

30 6. The method of claim 1, wherein the optical radiation is scanned from the source 

across the sample to simulate multiple sources. 
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7. The method of claim 1, wherein ten detectors are used to receive optical radiation 
from the sample. 



8. The method of claim 2, wherein the offset is from 0. 1 to 10 mm. 

9. The method of claim 3, wherein the one or more offsets are from 0.1 to 10 



mm. 



10. The method of claim 1, wherein the sample contains an absorbing object, and the 
image #/ a (r) is reconstructed by solving Equation 2a: 



A(r rf A) = -jjs Ma (r)L 0 (r„rM)G(T,n\ r d A d )d£ld 



V 



1 1. The method of claim 1 , wherein the sample contains a scattering object, and the 
image S/i s (r) is reconstructed by solving Equation 2b: 

A(r,A) = -JJ^,(r)[z^^ 

12. The method of claim 1, wherein the sample contains an absorbing and scattering 
object, and the image 8^{r) + ^ s(r ) is reconstructed by solving the summation of Equations 
2a + 2b. 



13. A system for indirect imaging of a sample comprising 

(a) a probe comprising a source optic fiber, one or more detector optic fibers, and a 
distal end, wherein a distal end of the source optic fiber and a distal end of each of the 
detector optic fibers extends through and ends in the distal end of the probe; 

(b) an optical radiation source connected with a proximal end of the source optic 

fiber; 

(c) one or more photodetectors, each connected to a proximal end of one of the 
detector optic fibers, to receive and convert optical radiation from each detector optic fiber 
into a digital signal corresponding to light emitted from the sample; and 

(d) a processor that processes the digital signal produced by the photodetectors to 
provide on an output device an image of spatially varying optical properties of one ■ 



: or more 
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features within the sample, wherein the image is reconstructed by performing a non-linear 
minimization between the digitized data and a transport-based photon migration model. 

14. The system of claim 13, wherein the distal end of each detector optic fiber is 
offset from the distal end of the source optic fiber. 

15. The system of claim 14, wherein the offset is from 0.1 to 10 mm. 

16. The system of claim 13, wherein the processor is programmed to process the 
digital signal to provide an image S Ma (r) that is reconstructed by calculating a matrix 
equation of the form y = Ax, wherein 

fit is the absorption coefficient, 
r is a position within the sample, 

5 X = 5fl a , 

y is the vector of detector signals for each source-detector pair, and 
Ay is an element of matrix A representing the product of (i) a Green function: 
G(r,n|r rf ,Q, ; ),foratransportequation(G): Mr) = /A + «0 and ^ s (r) = A*°. + 
S Ms (r), and (ii) the first-order (background) detector signal (Lo), for each measurement i at 
10 each voxel j. 

17. The system of claim 14, wherein one or more of the offsets of each detector- 
source pair are different. 

18. The system of claim 13, wherein the optical radiation is near-infrared radiation. 

19. The system of claim 13, wherein the optical radiation is scanned from the source 
optic fiber across the sample to simulate multiple sources. 

20. The system of claim 13, wherein ten detectors are used to receive optical 
radiation from the sample. 
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21. The system of claim 14, wherein the offset is from 0.1 to 10 mm. 

22. The system of claim 17, wherein the one or more offsets are from 0.1 to 10 mm. 

23. The method of claim 1, wherein the image is three-dimensional. 

24. The system of claim 13, wherein the image is three-dimensional. 

25. A probe for indirect mode imaging, comprising a source optic fiber, one or more 
detector optic fibers, and a distal end, wherein a distal end of the source optic fiber and each 
of the detector optic fibers extends through and ends in the distal end of the probe, and 
wherein the distal end of the source optic fiber is offset from each distal end of the detector 
optic fibers by 0.1 to 10 mm. 

26. The probe of claim 25, further comprising a scanning mirror to scan optical 
radiation across a sample. 

27. The probe of claim 25, wherein one or more of the offsets of one or more source- 
detector pairs are different. 
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